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Abstract 

Relativistic and non-relativistic ratios of Laplace transform QCD moment sum rules 
for charmonium are used in order to determine the value of the on-shell charm- 
quark mass. The validity of the non-relativistic version of QCD sum rules in this 
particular application is discussed. After using current values of the perturbative 
and non-perturbative QCD parameters, as well as experimental data on the J/ip 
system, we obtain rndQ'^ = rn^) = 1.46 it 0.07 GeV. 
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Recently, the value of the (on-shell) beauty-quark mass has been determined |I|] by con- 
fronting very accurate experimental data on the upsilon system with ratios of non- 
relativistic Laplace transform QCD moments. This theoretical framework, suggested by 
Bertlmann offers several advantages, e.g. radiative and non-perturbative corrections 
are well under control, and the non-relativistic limit follows quite naturally from quantum 
mechanical analogues [Q. This version of QCD sum rules leads to an expansion in powers 
of the inverse of the heavy quark mass which allows one to test the range of validity of 
the non-relativistic limit, and more generally, to assess the role of mass corrections. This 
might be of interest for calculations based on the simplifying assumption Kqcd/^q ^ 1- 

Non-relativistic Laplace moments appear to have a sensitive dependence on the quark 
mass. In fact, in spite of the large uncertainties affecting the values of Aqcd and the 
non-perturbative gluon condensate, mf, can be extracted from the upsilon data with high 
precision. This extraction is performed by confronting the ratios of Laplace transform 
moments calculated from experiment with those from theory. The latter involve the QCD 
parameters rrib, A, < a^G^ >, etc.. These ratios are functions of the Laplace variable, 
which acts as a short distance expansion parameter, and one finds a reasonably wide 
region in this variable where there is a matching between experiment and theory for a 
specific value of the quark mass. 

As pointed out in [|I|, a straightforward extension of this technique to the charm-quark 
may not work, as radiative and mass corrections could exceed 100 % . This would be 
true if the window in the Laplace variable would be the same for beauty and for charm. 
However, there is no a-priori reason for this to be the case. In fact, as also suggested in 
the matching between theory and experiment for the beauty and the charm quarks could 
take place at different ranges of the Laplace variable. If this range is such that radiative, 
non-perturbative, and mass corrections remain small, then it would become possible to 
extract the value of the charm quark from this framework. 
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In this note we study the ratios of relativistic Laplace transform QCD moments for the 
case of charm, and compare them with the non-relativistic versions. This provides a 
measure of the vahdity of the heavy quark mass expansion in this particular application, 
as well as an estimate of the systematic uncertainties affecting this technique. We make 
use of the current values of A and < agG"^ >, and carry out the non-relativistic expansion 
of the Laplace ratios to next-to-next to leading order in l/rric. We find this expansion to 
converge reasonably fast, and the radiative and non-perturbative contributions to be safely 
under control. Thanks to this feature and to the high accuracy of the experimental data, 
as in the case of m^, the extracted value of the (on-shell) charm-quark mass is affected 
by a relatively small uncertainty, in spite of the large uncertainties in A and < a^G^ >. 
In connection with the latter, we recall that most of the early (relativistic) QCD sum 
rule analyses of charmonium attempted to extract the values of both the charm-quark 
mass and the gluon condensate. Since < a^G^ > is now known independently from e'^e" 
and from r decay data |^, one is in a better position to determine m^. In fact, for a 
given pair of values of A and < a^G^ >, the matching between theory and experiment 
becomes a one parameter fit, the parameter being rric. We find that the values of rric 
from the fully relativistic ratios are in good agreement (within errors) with those from the 
non-relativistic ratios. We conclude with a comparison of our results with those obtained 
previously by other authors. 

We begin by considering the two-point function 

n^.(g) = t I rf^xexp(^gx)(0|T(V;(x)V;+(0))|0) = {-g,,q^ + g^g.)n(g2) , (1) 



with V^{x) = c(x)7^c(x). The function n(g^) has been calculated in perturbative QCD 
at the two-loop level H, with its imaginary part given by 



TT {v + 3) / IT 3 



2v A \2 Att 



[s-4m'J, (2) 



where v = Jl — Am'^/s, and is the charm-quark on-shell mass: rric = i^dQ^ = f^t] 
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The leading non-perturbative term in the operator product expansion of n(g^) involves 
the gluon condensate, i.e. 
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The function n(g^) satisfies a once-subtracted dispersion relation, and the subtraction 
constant can be disposed of by taking the Laplace transform 
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The quantity of interest to us here is the ratio of the first two Laplace moments, which 
can be expressed as 

d 

(5) 



n{a) = -—In n(a) 
acr 



From (2) and (3) one obtains 



n(cr) = exp(-4m^(T)7r A((t)[1 + a{a) + b{a) 



(6) 



where, with u = Am^a, 
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The function G{b,c,Lj) is related to the Whittaker function W\^fj,{uj) through 



G{b,c,uj)=u''-'/^e^/^Wx,^{u) 



(14) 



with /i 
result I 

where 



(c — 6) /2, and A = (1 — c — 6)/2. The ratio (5) can then be calculated, with the 
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The above expressions (15)-(18) involve no approximations, other than the two-loop per- 
turbative expansion, and retaining the leading non-perturbative term in the operator 
product expansion. We shall refer to Eq.(15) as the fully relativistic Laplace ratio. 

In the non-relativistic (heavy quark-mass) limit, the Laplace transform (4) becomes 

/■oo 

U(t) = / dEeM-TE) Im U(E) , (19) 
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where r = imccr, and s = (2mc + E^ so that E > 0. The Laplace ratio (5) is now given 

by 

(20) 



7^(T) = 2mc - — In U{t) . 
dr 

After expanding the functions G{b, c, a;) entering Eqs.(7)-(9), and (16)-(18), we obtain the 
non-relativistic ratio 
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(21) 



The appearance of ■s/rri^ above is only an artifact of the change of variables; written in 
terms of a, Eq.(21) contains no such term. The theoretical ratios of the first two Laplace 
moments (15) and (21) must now be confronted with a corresponding ratio involving 
the experimental data on the J/ip system. We parametrize the latter by a sum of two 



narrow resonances below DD threshold, followed by a hadronic continuum modelled by 
perturbative QCD, which gives 

27 1 1 r'^ 

n(a)Uxp = —J2T<y'Mvexp{-aMl.) + - / ds exp{-as) Im II{s)\qcd • (22) 

The experimental ratio is then calculated using (22) in (5). The continuum threshold 
So is chosen at or below the DD threshold. Reasonable changes in the value of sq have 
essentially no impact on the results, as n(cr) is saturated almost entirely by the first two 
J/ip narrow resonances. 

In the theoretical ratios we use the current values: A = 200 — 300 MeV , for four flavours 
0, and < ttsG"^ >= 0.063 — 0.19 GeV^ 0. We find that theoretical and experimental ra- 
tios match in the wide sum rule window: a ~ 0.8-1.5 GeV"^ for = 1.39-1.46 GeV 
in the fully relativistic case; and a ~ 0.6 - 0.8 GeV"^ = 1.40 - 1.53 GeV in the 
non- relativistic case. Figure 1 shows the behaviour of (15) for A = 200 MeV, and 
< asG^ >= 0.063 GeV^ (solid curve), together with the experimental ratio calculated 
from (22) (broken curve), corresponding to rric = 1.44 GeV. A similar qualitative be- 
haviour is obtained for other values of A and the gluon condensate, both in the rela- 
tivistic and the non-relativistic versions of the sum rules. For values of cr inside the sum 
rule window, the hierarchy of the various terms in the non-relativistic Laplace ratio (21) 
guarantees a fast convergence. In fact, the leading correction in l/rric is at the 15-20% 
level, the radiative correction and the non-perturbative contribution amount both to less 
than 10% . At the same time, the next, and next-to-next to leading (in l/rric) terms 
everywhere in (21) are safely small, as it can be easily verified from (21) noticing that if 
o" ^ 1/2 GeV~^, then r ~ 2mc. Clearly, the complete analysis at the level of accuracy of 
these next-to-leading mass corrections would require the evaluation of the perturbative 
0{a1) terms. Combining the results from both versions of the Laplace ratios, leads to the 
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result 



m^(Q2 = = ± 0.07 GeV 



(23) 



In order to facilitate the comparison of (23) with previous determinations based on various 
versions of QCD sum rules 0) (H - IH? "we show in Fig. (2) the dependence of nic on A 
for three different values of < a^G^ >, and in Fig. (3) the dependence of rric on the gluon 
condensate for A in the range: A = 100 — 400 MeV. Both figures correspond to the fully 
relativistic version of the QCD sum rules. Figures 4 and 5 refer to the non-relativistic 
determination. In comparing values of rric from different determinations, it is important 
to know which values of A and < a^G^ > have been used, as well as which renormalization 
point has been chosen, e.g. some authors determine rndQ'^ = —ml), which is related to 
the on-shell mass mc(m^) through 

mliml) = ml{~ml){l + ^^a,) . (24) 

TT 



After using the same values of A and < asC^ > as used in ||^, and after converting 

to the on-shell mass (if necessary), we find that results from the present method are in 
very good agreement with those of and [Q, and agree within errors with the 
latter determinations being on the low side of our error bars. On the other hand, the 



technique used here gives values of somewhat higher than those obtained in |ll2|-|14 



Recently [jT3[, the two- loop correction to the Wilson coefficient of the gluon condensate 
has been calculated. We have incorporated this additional term in the Laplace ratios, 
and find that the correction it introduces is about a factor of 2 larger than the l/rric 
correction to < a^G^ > but with an opposite sign. Given the relative smallness of the 
overall contribution from the gluon condensate, and the conservatively large uncertainty 
we have allowed in its value, the final result for rric remains basically unchanged. 

Finally, regarding the recent observation that the definition of heavy quark pole mass in 
the context of the Heavy Quark Effective Theory (HQET) should contain some intrinsic 
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ambiguity beyond perturbation theory [|16] , due to nonperturbative long distance effects, 
we notice that the hadronic system considered here is made of two heavy quarks, while 
HQET strictly applies to heavy-light bound states.^ The non-relativistic ratio (21) is not 
a relation of the HQET, although it is obtained formally from the relativistic ratio (15) 
in the large quark-mass limit. The purpose of considering (21) together with (15) has 
been to assess the size of relativistic corrections, as well as of systematic uncertainties of 
QCD sum rules. The numerical results, and in particular the consistency between the 
relativistic and nonrelativistic determinations of 171^, indicate that these effects should be 
small. 
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Figure Captions 

Figure 1: The fully rclativistic ratio (15) (solid curve), and the experimental ratio (broken 
curve). The values of the parameters are: A = 200 MeV, {as G^) = 0.063 GeV^, 
and rric = 1.44 GeV. 
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Figure 2: Dependence of rUc on A for (a, G^) = 0.038 GeV^ (curve (a)), 0.063 GeV^ (curve 
(b)), and 0.19 GeV"^ (curve (c)). The fully relativistic ratio (15) has been used. 
Curve (a) is shown for reference purposes, as this low value of {ag G^) was not used 
in our analysis. 

Figure 3: Dependence of rric on {ag G'^) for A = 100 MeV (curve (a)), 200 MeV (curve (b)), 
300 MeV (curve (c)), and 400 MeV (curve (d)). The fully relativistic ratio (15) has 
been used. Curves (a) and (d) are shown for reference purposes, as these extreme 
values of A were not used in our analysis. 

Figure 4: Same as in Fig. 2, except that the non-relativistic ratio (21) has been used. 

Figure 5: Same as in Fig. 3, except that the non-relativistic ratio (21) has been used. 
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